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Verfahren zum Orten von Lecksteilen in Rohrleitungen 

Bei dem erfmdungsgema&en Leckortungsverfahren mit 
Kreuzkorrelationsanalyse wird zur Messung der Laufzeitdtf- 
ferenz ein 2-Kanal-FFT-Analysator benutzt. Das Kreuzlei- 
stungsspektrum wird zunachst einer Egalisierung unterwor- 
fen, erforderlichenfafls mit einer geeigneten Frequenzbe- 
wertung gewichtet und anschliefcend in den Zeitbereich zu- 
rucktransformiert. Dadurch wird an Stella der bekannten 
Kreuzkorreiatjonsfunktion eine andere Zeitfunktion mit op- 
timalem Auflosungsvermogen erhajten, was die Ermittlung 
von Lecksteilen mit hdherer Genauigkeit und in ein em gro- 
Ceren MeRbereich ermoglicht. Die nachtragliche Egalisie- 
rung des Kreuzleistungsspektrums ist in viel feinerer Weise 
moglich, als sie mit konventionellen Filtern 3 nach Figur 1 
erreichbarware. 

Das Verfahren laBt sich ohne zusatzlichen Gerateaufwand in 
ublichen FFT-Analysatoren implementieren. Da es automa- 
tisch ablaufen kann, erfordert die Zwischenschaltung der 
Egalisierung kein en zusatzlichen Zeitauf wand im prakti- 
schen Einsatz. 
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Patentanspruche kann dann aus der gemessenen Laufzeitdifferenz r m der 

Signale a(t) und b(t) der Ort der Leckstelle berechnet 

1. Verfahren zum Orten von Leckstellen in Rohrlei- werden. 

tungen, insbesondere in Wasserversorgungsnetzen, Als MeBgerat 4 fur die Laufzeitdifferenz wurden zu- 

durch Messung der Laufzeitdifferenz der von ei- 5 nachst Korrelatoren verwendet, die im Zeitbereich ar- 

nem Leek ausgehendeh Schallwellen, wobei von beiten und die Kreuzkorrelationsfunktion 

akustischen Sensoren gelieferte elektrische Signale ______ 

durch Kreuzkorrelation mittels eines Zweikanal- RBA(x)*= a(t) • b(t+x) 
FFT-Analysators ausgewertet werden, dadurch ge- 

kennzeichnet, daB das im Analysator gemessene 10 liefern. Beschrieben ist eine solche MeBanordnung in 

komplexe Kreuzleistungsspektrum zunachst in sei- Spalte 6 Zeilen 14 bis 23 und 37 bis 49 der DE-PS 

nem Betrag durch Abschwachung von Frequenz- 30 45 660. 

bereichen mit groBer Amplitude und Anheben von Neuerdings sind aber auch FFT-Analysatoren gemaB 

Frequenzbereichen mit geringer Amplitude so ge- Fig. 2 im Einsatz, bei denen zunachst die Zeit-Signale 

formt wird, daB es in einem moglichst breiten Fre- 15 a(t) und b(t) mit der Fourier-Transformation in komple- 

quenzband moglichst gleiche Amplituden aufweist xe Spektren A(w) und B(<o) uberfiihrt werden: 
und erst nach dieser Egalisierung die Rucktransf or- 

mation in den Zeitbereich vorgenommen wird. A(w)=* F{a( rj}, B(o>)= F{b( fj}. 
Z Verfahren nach Anspruch 1, dadurch gekenn- 

zeichnet, daB die Egalisierung durch Dividieren des 20 AnschlieBend wird das Kreuzleistungsspektrum 

Kreuzleistungsspektrums mit der Wurzel aus dem _________ 

Produkt der Leistungsspektren der zwei Eingangs- GBA(w)= A(o>/ • B(<oJ 
signale erfolgt 

3. Verfahren nach Anspruch 2, dadurch gekenn- (A(q>)* ist die konjugiert komplexe GroBe zu A(o)J) be- 
zeichnet, daB im egalisierten Kreuzleistungsspek- 25 rechnet und danach die Kreuzkorrelationsfunktion 
trum eine zusatzliche Frequenzbewertung, vor- durch die Invers-Fourier-Transformationgewonnen: 
zugsweise in Form einer vom Benutzer variierba- 

ren Filter-Charakteristik, vorgesehen wird, um RBA(r)= F- l {GBA(&$. 
eventuelle Storanteile, welch e bedingt durch die 

Egalisierung eine groBere Amplitude als im ur- 30 Bei c^er Laufzeitdifferenz x m erscheint in RBA ein Ma- 

sprunglichen Kreuzleistungsspektrum aufweisen ximum. Je scharfer dieses Maximum ausgepragt ist, de- 

konnen, vor der Rucktransformation in den Zeitbe- sto sicherer und genauer kann r m bestimmt werden. Ein 

reich zu eliminieren. rasches Abklingen von RBA um r m herum ist besonders 

auch dann wichtig, wenn mehrere Schallquelten mit ver- 

Beschreiburig 35 schiedenen Laufzeitdifferenzen vorhanden sind. 

Wie rasch RBA um r m herum abklingt, wird durch die 

Die Erfindung betrifft ein Verfahren zum Orten von Hullkurve von GBA bestimmt: schreibt man das kom- 

Leckstellen in Rohrleitungen, insbesondere in Wasser- plexe Kreuzleistungsspektrum nach Betrag und Phase, 

versorgungsnetzen, durch Messung der Laufzeitdiffe- so folgt fur RBA: 

renz der von einem Leek ausgehenden Schallwellen, 40 

wobei von akustischen Sensoren gelieferte, elektrische RBA(x)*= F-*{\GBA(<d)\}* F-*{exp(j0%. 
Signale durch Kreuzkorrelationsanalyse ausgewertet 

werden. Eine MeBanordnung fur ein solches Verfahren Der Stern bezeichnet die Faltungsoperauon [2], Seite 

zeigt die Fig. 1. An den Orten a und b sind im Abstahd / 16. Die gesuchte Information iiber die Laufzeitdifferenz 

an der Rohrleitung die akustischen Sensoren 1 ange- 45 r„, ist ausschlieBlich in der Phase 0 enthalten. Es ist 

bracht, deren elektrische Signale mit Verstarkern 2 und daher zulassig, auf den Term links vomFaltungssymbol, 

Filtern 3 aufbereitet werden. Mit den Filtern 3 werden der Iediglich die abklingende Schwingung von RBA um 

storende Signalanteile, z. B. das tieffrequente Pumpen- r m herum beschreibt, EinfluB zu nehmen. Ziel dieser 

gerausch, vom Analysator 4 ferngehalten. Zur Optimie- Operation ist es, eine ausgeglichene und breitbandige 

rung der Filtereinstellung Iiefert nach [1] die KohaYenz- 50 HQllkurve des Kreuzleistungsspektrums zu erzieien, 

funktion wertvolle Zusatzinformation. Nach [2], Seite weil die Auflosung im Zeitbereich At gemaB der Un- 

54 ist scharferelation B • At=\ der Bandbreite B der zugeho- 

rigen Frequenzfunktion umgekehrt proportional ist Die 

y 2 *=\GBA\ 2 /(GAA • GBB) % Notwendigkeit f Qr eine derartige Verbesserung der bis- 

55 herigen MeB verfahren wird in Fig. 3 mit MeBergebnis- 

wobei mit GBA das Kreuzleistungsspektrum, mit GAA sen bei der Leckortung auf einer Fernwarmeleitung 

und GBB das Leistungsspektrum des Sensorsignals in (Stahlrohr, Nennweite 100 mm, Lange des MeBab- 

Kanal a und b bezeichnet wird. Danach werden die Si- schnittes 195 m) demonstriert: In den Spektren der Sen- 

gnaie dem Analysator 4 zur Bestimmung der Laufzeit- sorsignale erscheinen scharfe Resonanzen (wenig be- 

differenz zugefuhrt. 60 dampfte Biegeschwingungen des Stahlrohrs), die sich in 

Unterhalb einer durch den Durchmesser des Rohres das Kreuzleistungsspektrum abbilden. Daher erscheint 

gegebenen Grenzfrequenz sind im Fluid nur ebene in RBA nur ein breites verschmiertes Maximum, die 

Schallwellen ausbreitungsfahig. Unter dieser Vorausset- Bestimmung von r m ist mit groBer Unsicherheit behaf- 

zung ist die Laufzeit, die das vom~Ausstr6mvorgang im tet Zur Vermeidung dieses Nachteils wurde-bisher nur 

Leek L hervorgerufene Gerausch benotigt, um zu den 65 versiicht, durch eine geeignete spektrale Formung in 

Sensoren a oder b zu gelangen, der Entfernung direkt den Filtern 3 ein ebenso ausgeglichenes Spektrum der 

proportional. Bei bekannter Schallgeschwindigkeit — Sensorsignale vor dem KorrelationsprozeB herbeizu- 

sie wird in der Regel in einem Vorversuch bestimmt — fiihren und so die Korrelationsfunktion zu verbessern 
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(United States Patent 44 35 974)._Diese Vorgehensweise 
verlangt zusatzlichen Cerate- und Zeitaufwand und 
fuhrt in der Anordnung nach Fig. 1 zu einem zusatzli- 
chen MeBfehler, weil sich bei Einstellung auf zwei ii 
verschiedene Spektren A(a>) und B(o>) die Filter in ihrem 
Phasengang unterscheiden. 

Der Erfindung liegt deshalb die Aufgabe zugrunde, 
dieses aus [1] bekannte, mit einem FFT-Analysator ar- 
beitende Leckortungsverfahren- so zu verbessern, dali 
ohne zusatzlichen Gerate- und Zeitaufwand eine ge- 
nauereLeckortungerm6gIichtwird _ 

Die Losung dieser Aufgabe geschieht durch die im 
Kennzeichen des Anspruchs 1 angegebenen Verfan- 
rensschritte. Im folgenden wird die Erfindung anhand 
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der Figuren erlautert Es zeigt 

f ig. i die aus ucr ul^-i ^ ^ -w v/ - ,v ~~* v — •— 

nU Rg. 2 den anstelle des MeBgerates 4 in Fig. 1 ver- 
wendbaren FFT-Analysator, 

Fig. 3 ein MeBbeispiel an einem Rohrabschmtt mit 20 
dembisherublichenFFT-Verfahren. 

Fig. 4 eine Gegenflberstellung von MeBergebnissen 
nach dem bisherigen und dem erfindungsgemaBen FFT- 

Verfahren, . , , , -~ 

Fig. 5 Ergebnisse mit dem Verfahren nach der Erfin- 
dung an dem der Fig. 3 zugrunde liegenden Ronrab- 

^^e^angestrebte Egalisation des Kreuzleistungsspek- 
trums GBA wird erfmdungsgemaB rhit der sogenannten 
DissipationsfunktionZ>Fdurch No Wrung auf die JLei- 
stungsspektren der Eingangssignale GAA und GBB 
nach folgender Definitionsgleichung erreicht: 
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DF(o>) 



GBA(g>) 



VGAA(o)'GBB(o)J 



- yexp(/0) 



Laufzeitdifferenzen (Leek nahe der Mitte des MeBab- 
schnittes)gemaB^f«l/Berheblicherweitert. 

Bedingt durch die endliche Mittelungszeit bei der 
Messung des Kreuzleistungsspektrums und durch das 
Auftreten instationarer akustischer Storsignale ver- 
schwindet GBA auch auBerhalb des Koharenzbereiches 
nicht vollig, sondern weist dort einen zumeist regellosen 
Verlauf mit sehr kleinen Amplituden auf. Durch die 
Normierung von GBA auf die Spektren der Eingangssi- 
gnale kann der Fall eintreten, daB solche Bereiche in 
GBA stark hervorgehoben werden. und in der Dissipa- 
tionsfunktion als Storanteil erscheinen. Fig. 5 illustnert 
diesen Vorgang fur das Beispiel in Fig. 3. In diesem Fall 
ist es sinnvoll, durch eine zusatzliche Gewicntung von 
PFdiese Storanteile vor der RUcktransformaUon in den 
7eitbereich auszublenden. Fig. 5 zeigt die bewertete 
DF W wobei fur den Obergang vom DurchiaJioereicn 
zum Sperrbereich in Anlehnung an das bekannte Han- 
ning-Fenster [2\ Seite 74, cos 2 -Flamcen verwendet wur- 
den. Im Vergleich zu GBA in Fig. 2 ist ersichthch, daB in 
der DF jetzt mehrere Oszillationen etwa gleicner Am- 
plitude vorhanden sind, die aus DF W gewonnene Dissi- 
pationszeitfunktion DTF W weist ein zweifelsfreies Maxi- 
mum auf. 

ZitierteLiteratur • 
[1] Fuchs, H.V.; Riehle. R.; Schmidt-Schykowski, 
KL, Systematische Rohrnetz-Oberprufung mit Hilfe 
verbesserter akustischer Methoden in: Fortschritte 
der Akustik - DAGA '85, Seiten 739-742, Bad 
Honnef:DPG-GmbH1985. . 
[2] Bendat, J.S.; Piersol, A.G., Engineering Appli- 
cations of Correlation and Spectral Analysis, John 
Wiley & Sons, New York 1 980. 
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y2 bezeichnet wieder die Koharenzfunktion der Si- 
gnale aft) und b(t). Die Dissipationszeitfunktion DTF 
wird durch die Invers-Fourier-Transformation aus der 
Dissipationsfunktion gewonnen: 

DTF(t)= F-i{DF(co)i - F-^yJ* F-i{exp(j0fr 

Solange die Sensorsignale iiberwiegend durch das 
Leckgerausch bestimmt sind, liegt der Wert der Koha- 45 
renzfunktion nahe bei 1. In diesem Bereich 1st y prak- 
tisch unabhangig vom Leistungspektrum der Signal- 
quelle und vom Ampjitudengang des Obertragungswe- 
ges vom Leek zu den Sensoren. Daher weist die Koha- 
renzfunktion fast immer eine groBere Bandbreite auf als 50 
das Kreuzleistungsspektrum - und demgemaB klingt 
die Dissipationszeitfunktion urn v m herum schneller ab 
als die Kreuzkorrelationsfunktion. Fig. 4 veranschau- 
iicht diese Verhaltnisse an Hand von MeBergebnissen 
fureinl2mlangeswassergefailtesStahlrohr 0 1 .Ge- 55 
rauschquelle ist der Strdmungsvorgang durch eine 
Lochblende. Er erzeugt ein typisches Leckgerausch mit 
uber der Frequenz abfallendem Spektrum. Auch die 
Hullkurve des Kreuzleistungsspektrums 1st von dieser 
Form. GBA ist nur in einem schmalen Band bei tiefen eo 
Frequenzen ausgepragt, daher weist RBA urn die ( Werte 
-r m und + r,„ (einlaufende und reflektierte Welle) eine 
relativ langsam abklingende Oszillation tiefer Frequenz 
auf. Demgegenuber ist die DissipatLonsfunktion DF in 
einem sehr breiten Frequenzband gemaB der Koha- 65 
renzfunktion ausgepragt, so daB die Dissipationszeit- 
funktion ein sehr scharfes Maximum bei -r iW ergibt. 
* Zugleich wird der MeBbereich in Richtung sehr kleiner 
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